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Abstract 

Let A be a max-stable random vector with positive continuous density. It is 
proved that the conditional independence of any collection of disjoint subvectors of X 
given the remaining components implies their joint independence. We conclude that 
a broad class of tractable max-stable models cannot exhibit an interesting Markov 
structure. 
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1 Introduction 

As pointed out by Dawid (1979) independence and conditional independence are key con¬ 
cepts in the theory of probability and statistical inference. A collection of (not necessarily 
real-valued) random variables Yi,...,!^ on some probability space (tl,^,P) are called 
conditionally independent given the random variable Z (on the same probability space) 
if 

k 

P(yi e Ai, ..., Yfc G Afc I Z) = [] P(yi €Ai\Z) P-a.s., 

i=l 

for any measurable sets Ai,..., from the respective state spaces. The conditioning is 
meant with respect to the cr-algebra generated by Z. A particularly important example for 
the conditional independence to be an omnipresent attribute are the Gaussian Markov 
random fields that have evolved as a useful tool in spatial statistics (Lauritzen 1996, 
Rue &: Held 2005). Here, the zeroes of the precision matrix (the inverse of the covariance 
matrix) of a Gaussian random vector represent precisely the conditional independence of 
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the respective components conditioned on the remaing components of the random vector. 
Hence, sparse precision matrices are desirable for statistical inference. 

In the analysis of the extreme values of a distribution (rather than fluctuations 
around mean values) max-stable models have been frequently considered. We refer to 
Blanchet Sz Davison (2011), Buishand et al. (2008), Engelke et al. (2014), Naveau et al. 
(2009) for some spatial applications among many others. Their popularity originates 
from the fact that max-stable distributions arise precisely as possible limits of location- 
scale normalizations of i.i.d. random elements. A random vector X is called max-stable 
if it satisfies the distributional equality anX -\-bn = max(A(^\..., for indepen¬ 

dent copies X^^\ ..., of X for some appropriate normalizing sequences > 0 and 
bn € K, where all operations are meant componentwise. If the components Xi of X are 
standard Frechet distributed, i.e. P(Aj < x) = exp(—1/x) for x G (0, oo), we have an = n 
and bn = 0 and the random vector X will be called simple max-stable. 

Let / be a non-empty finite set. It is well-known (cf. e.g. Resnick (2008)) that the 
distribution functions G of simple max-stable random vectors X = (Aj)jg/ are in a one- 
to-one correspondence with Radon measures H on some reference sphere = {cu € 
[ 0 ,oo)^ : ||a;|| = 1} that satisfy the moment conditions f a;jLf(dw) = 1, f G /. The 
correspondence between G and H is given by the relation 

G(x) = F(Xi < Xi, i € I) = exp ( — max — Lf(dti;) ) , xG( 0, oo)^. 

V Js+ Xi ) 

Here, H-H can be any norm on and H is often called angular or spectral measure. 

In general, neither does independence imply conditional independence nor does condi¬ 
tional independence imply independence of the subvectors of a random vector. Consider 
the following two simple examples which illustrate this fact in the case of Gaussian ran¬ 
dom vectors (Example 1) and max-stable random vectors (Example 2). For notational 
convenience, we write X ALY if X and Y are independent and X ALY \ Z if X and Y 
are conditionally independent given Z and likewise use the instructive notation -LL(T^ Xi 
and Xi\ Z if more than two random elements are involved. 

Example 1. Let X\,X 2 ,X^ be three independent standard normal random variables 
and, moreover, X 4 = Xi -\- X 2 and X^ = Xi A- X 2 A- X^. Then all subvectors of (Aj)f^^ 
are Gaussian and 

Xi _LL X 2 , but not Xi _LL X 2 \ X^, (1) 

whereas Xi _LL X 5 | A 4 , but not Xi _LL A 5 . (2) 

Example 2. Let Xi,X 2 ,X^ be three independent standard Frechet random variables 
and, moreover, X 4 = msix{Xi, X 2 ) and A 5 = max(Ai, Ai 2 , A 3 ). Then all subvectors of 
(Aj)j=i _...^5 are max-stable and both relations ( 1 ) and ( 2 ) hold true also in this setting. 

However, if the distribution of a max-stable random vector has a positive continuous 
density, then conditional independence of any two subvectors conditioned on the remain¬ 
ing components implies already their independence. To be precise, when we say that a 
random vector has a positive continuous density, we mean that the joint distribution of its 
components has a positive continuous density. The following theorem is the main result 
of the present article. If A = (Aj)jg 7 is a random vector, we write Xa for the subvector 
{Xi)i^A if A C I. The same convention applies to non-random vectors x = {xi)i£j. 
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Theorem 1. Let X = (Xj)jg/ he a simple max-stable random vector with positive con¬ 
tinuous density. Then, for any disjoint non-empty subsets A and B of I, the conditional 
independence Xa AL Xb \ Xj\(^aub) implies the independence Xa AL Xb- 

A proof of this theorem is given in Section 3. Beforehand, some comments are in 
order. 

(a) First, the requirement of a positive continuous density for X is much less restrictive 
than requiring the spectral measure Ff of X to admit such a density, cf. Beirlant et al. 
(2004) pp. 262-264 and references therein. For instance, fully independent variables X = 
{Xi)i^j have a discrete spectral measure, while their density exists and is positive and 
continuous. A more subtle example is, for instance, the asymmetric logistic model (Tawn 
1990), which admits a continuous positive density and whose spectral measure carries 
mass on all faces of 5+, cf. also Example 3. 

(b) Secondly, both random vectors (Xj)j=i_ 2,5 and {Xi)i=iA,b that were considered in 
the Gaussian case in Example 1 have a positive continuous density on Hence, there 
exists no version for Theorem 1 for the Gaussian case. 

(c) Note that the implication of Theorem 1 is the independence of Xa and Xb, not 
the independence of all three subvectors Xa-,Xb,Xj\^^avjb)- 

(d) By means of the same argument that shows that pairwise independence of the 
components of a max-stable random vector implies already their joint independence, we 
may deduce a version of Theorem 1, in which more than two subvectors are considered. 

Corollary 2. Let X = [Xi)i£i he a simple max-stable random vector with positive contin¬ 
uous density. Then, for any disjoint non-empty subsets Ai,... ,Ak of I, the conditional 
independence | ^ implies the independence 

(e) The non-degenerate univariate max-stable laws are classified up to location and 
scale by the one parameter family of extreme value distributions indexed by 7 G M 


F-y{x) = exp(-(l -I- jx) 


X G < 


(-1/7,00) 


(-00,-1/7) 


7 > 0, 
7 = 0, 
7 < 0. 


Any other (not necessarily simple) max-stable random vector is obtained through a trans¬ 
formation of the marginals that is differentiable and strictly monotone on the respective 
sub-domain on (cf. e.g. Resnick (2008) Prop. 5.10). Hence, the above results remain 
valid for the general class of max-stable random vectors. 

(f) Dombry &: Eyi-Minko (2014) show that, up to time reversal, only max-auto¬ 
regressive processes of order one can appear as discrete time stationary max-stable pro¬ 
cesses that satisfy the first order Markov property. This result indicates already that 
the conditional independence assumption is to some extent unnatural in presence of the 
max-stability property. 


Example 3. Various classes of tractable max-stable distributions admit a positive con¬ 
tinuous density, such that Theorem 1 and Corollary 2 apply. Popular models that are 
commonly used for statistical inference include the asymmetric logistic model (Tawn 
1990), the asymmetric Dirichlet model (Coles &: Tawn 1991), the pairwise beta model 


3 



(Cooley et al. 2010) and its generalizations involving continuous spectral densities (Ballani &: Sclilather 
2011) in the multivariate case. Moreover, most marginal distributions of spatial models 
such as the Gaussian max-stable model (Genton et al. 2011, Smith 1990) or the Brown- 
Resnick model (Hiisler & Reiss 1989, Kabluchko et al. 2009) possess a positive continuous 
density if the parameters are non-degenerate. Hence, if any of the components of the pre¬ 
viously mentioned extreme value models exhibit conditional independence given any of 
the remaining components, they must be independent. 

In the remaining article we subsume auxiliary arguments in Section 2 and give all 
proofs in Section 3. 

2 Preparatory results on max-stable random vectors 

Throughout this section let G be the distribution function of a simple max-stable random 
vector X = {Xi)i^i that has a positive continuous density. We denote its exponent 
function by 

V(x) := —logG(x) = [ max ( — ) a:e(0, cx))^. 

Js+ \xij 

Lower order marginals that refer to a subset A of I are obtained as minjg^c(xi) —>■ oo, 
where = I \ A. We write —>■ oo for —>■ oo, and with this notation 

G^{xa) ■= lim G{x) and V^{xa) ■.=—\ogG^{xA)- 

Since G is absolutely continuous, the partial derivatives 

B\b\ A\b\ 

Gs(xa) ■= -^G^{xa) and Vb{xa) ■= 

exist and are continuous for B C A, and the latter are homogeneous of order —(|il|-|-l) 

(Coles & Tawn 1991). An elementary computation shows that 

Gb{xa) = Wg(a;A)exp {-V^(xa)) , 


where 


wii{xM)= 

7rgn(M) JGtt 

and n(M) stands for the set of partitions of M for M <Z N d I. 

Let us further denote the set of non-empty subsets of I by C{I). The collection of 
exponent functions {V^)a£C{I) is in a one-to-one correspondence with its Mobius inversion 
{dA)AeC{i), i-e-, if we set 

dA{x) := 


4 




it follows that can be recovered via 

V^{xa) = ^ dsix) (3) 

(cf. Papastathopoulos &: Tawn (2014), Theorem 2 and Schlather & Tawn (2002), Theo¬ 
rem 4 or, more generally, Berge (1971) Chapter 3, Section 2 for the Mobius inversion). 
Finally, we define 

Xa{xa) ■■= lim dA{x) = V (-1)I^I+V^(xb) = dsix). 

X^C^OO ^^ ^^ 

BgC{ 1):B(ZA BGCiI):A(ZB 

Then the collection of functions {xa)agC{i) is also in a one-to-one correspondence with 
{y^)AeC(^i) well as {dA)AeC{i) the inversions are given by 

dA{x) = E {-'i^)\^\A\xB{xB), 

B£C{I):AcB 

V^{xa) = {-1)\^\"^^xb{xb)- 

BeCiiy.BcA 

Further expressions for V^, dA and XA are collected in Lemma 3. Note that xa{xa) > 
dA{x) and thus. 


= 0 = 0 . 


( 4 ) 


Lemma 3. The functions and dA and XA (with A G C{I)) can be expressed in terms 
of the spectral measure H as follows: 


V^{xa) = / max 


'S+ 


OJi 


dA{x) = 

Js+ 


i&A \ X 
LOi 

mm ■ 
i&A \Xi 


Xa{xa) = I min ( — ) H{duj). 

Jg+ y Xi 


H{duj) 


— max ■ 
jeA^ \Xj 


Here = max(0, z) and max(0) = 0. 


H{dLj), 


It turns out that the following two quantities are closely linked to conditional indepen¬ 
dence and independence of subvectors of X, respectively. For non-empty disjoint subsets 
A,B of I and C = I \ (^4 U B), we set for x G (0, oo)^ 

dA,B{x) := V^^'^{xA\jc) + V^^^{xB\jc) -V{x) - V^{xc) 

= 

LGC(7'):LnA7^0,LnS^0,LnC=0 

XA,b{xaub) ■= lim dA,B{x) = V^{xa) + V^{xb) - V^^^{xaub) 

XC^OO 

= 

LeCii):LnA^(li,LnBj^<l) 


where the last equalities follow from (3). The following lemma is an analogue of Lemma 3. 
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Lemma 4. The functions ci'^d xa,b can be expressed in terms of the spectral measure 
H as follows 


dA,B{x) = 

Xa,b{xaub) = 


Xi 


— max 


. f (Ui\ ( UJi 

iin max — , max — 

\i£A \Xi J i£B \Xi J J jeC 

[ min fmax f^ , max ^ ^ H{duj). 
Js+ \i&A \XiJ i&B \XiJ J 




Note that xa,b{xavjb) > dA,B{x) implies similarly to (4) that 


dA,B =0 XA,B = 0. 


(5) 


General expressions for the regular conditional distributions for the distribution of a 
max-stable process conditioned on a finite number of sites that are based on hitting scenar¬ 
ios of Poisson point process representations have been computed in Dombry &: Eyi-Minko 
(2013), Oesting (2015), Oesting & Schlather (2014) under mild regularity assumptions or 
in Wang & Stoev (2011) for spectrally discrete max-stable random vectors. 

Let again A and B be non-empty disjoint subsets of I. Since we assumed a pos¬ 
itive continuous density for G (and hence also for its marginals), the numerators and 
denominators in 


G{xa\xb) ■= 


G"^^^ (xaub) 
G^ixB) 


exp (- [E^^^(xaub) - y^{xB)]) 


W^^^{xaub) 

WE{xb) 


are non-zero and continuous for x G (0, oo)^ and the expression G{xa\xb) constitutes a 
regular version of the conditional probability ^{Xa < xa\Xb = xb)- 


Proposition 5. The functions XA,B and dA,B are connected with the independence and 
conditional independence of the respective subvectors of X as follows. 

a) Xa -LL Xb I Xj\(^aub) dA,B = 0. 

b) Xa -LL Xb ^ xa,b = 0. 

Remark. The assumption that G admits a positive continuous density on (0, oo)^ is crucial 
for part a) to hold true. It fails in Example 2. 

Moreover, it is a simple consequence of Berman (1961/1962) and de Haan (1978) that 
the pairwise independence of any disjoint subvectors of the simple max-stable random 
vector X implies already their joint independence. 


Lemma 6. If Xai, ■ ■ ■ ,Xa^, are pairwise independent subvectors of a simple max-stable 
random vector X (for necessarily disjoint Ai C //, then they are jointly independent. 
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3 Proofs 


Proof of Lemma 3. The first equation is clear from the definition of V^. The relation for 
cIa can be obtained as follows. 

dA(x) = 

BeC(I):A‘=CB 




B&C{I):A‘^CB 
UJi 

mm I — I — max 


a;,: 


H{duj) 


S+ \_i&A \Xi J j&A<^ \XjJj_^_ 


IXn 


H{duj) 


In order to obtain the last equality, we denote Oj = Ui/xi and distinguish two cases: 
1st case: A = I. Then 

(—max (oj) = (—max (uj) = min (oj) . 

^^ i£B ^^ zGB i£l 

B&C{I):A‘^(ZB B€C{1) 

2nd case: A ^ I. Then set b := maxjg^c Oj and Cj := max(ai, b), such that 


(-l)l®n^l+i max {a) - b 
ieB ' ^ ^ ieBnA^ 

BeC{l):A^(ZB BeC{l):A^CB,B^A<^ 


max (cj) — b = min (cj) — b 
i&J ic A 


UCA:U^0 


eA 


lin (maxioj, b)) — b = max ( min(oi), 6 ) — ^ = ( min(aj) — b ) 
eA \ieA J \i&A J 

The expression for XA follows immediately. 


= mm I 
ieA 


□ 


Proof of Lemma 4- Similar to the proof of Lemma 3, the relation for dA,B follows from 
dA,B{^) = V^^^ixAuc) + V^^^ixBuc) - Vix) - V^ixc) 


OJi 


UJi 


UJr 


UJi 


/ max ( — ) + max ( — 1 — max ( — ) — max ( — ) Hiduj) 
Js^i&AUC \XiJ i£BUC\XiJ \Xi/ - 


f -I I 

/ mm max — , 

Js+ I \XiJ 


i&B \ X. 


UJi 

max I — I I — max 


jec \x 


i&C \Xi_ 

H{duj), 


^3 / J 


where the last equality is obtained from 


max (a,) + max (a, ) — max (a,) — max (a,) 
iG^uc^ ^ ieBuc^ ’ iei ^ ’ ieiC ^ ’ 

= min I max (oj), max (oj) ) — max(ai) 


ieAuc 


ieBuc 


iec 


= max 


mm 


( inax (oj), max (oj) ) , inax (oj) 
V i(^A i^B J i^C 


— max (oj) 
i&C 


/ 



min ( max (a, ) 

, max (oj) 

— max (oj) 

\ieA^ 

i&B ' J 

idC 


if we denote ai = ujijxi. The expression for XA,B follows immediately. 


□ 
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Proof of Proposition 5. a) As before, let C = I\{A[JB). Since G{x) = exp(—1/(x)) has a 
positive continuous density, we have that the conditional independence Xa Pl Xb \ Xc 
for C = I \ (A U -B) implies that for all x G (0, ooY 

G(xa\Xc) G{xb\Xc) = G{xaub\Xc) P-a.s.. 


Since Xq has a positive continuous density with respect to the Lebesgue-measure on 
(0,oo)^, it follows that 


G{xa\xc) G{xb\xc) = G{xa\jb\xc) for all x G Q, 


where Q is a dense subset of (0,oo)^. By the continuity of these expressions in x G 
(0, coY, the equality holds for all x G (0, coY and is equivalent to 


exp {dA,B{x)) 


WY^(^{xa,jc)WS^^{xbuc) 

WY^b^(^{xavjbuc)WS{xcY 


X G (0, ooY ■ 


( 6 ) 


Here, (Ia^b > 0 and dA,B is homogeneous of order —1, while the components Vj^ that 
build the terms are homogeneous of order —(| J| + 1). Now, replacing x by t~^x 
for t > 0 in (6), we see that the left-hand side grows exponentially in the variable t 
as t tends to oo if dA,B{x) > 0, while the right-hand side exhibits at most polynomial 
growth. Therefore, dA,B{x) = 0 for x G (0, oo)^. 


b) Both sides are equivalent to G^{xa)G^{xb) = G^^^{xaub) for all x G (0,oo)-^. □ 

Proof of Theorem 1. The hypothesis follows from Proposition 5 and (5). □ 


Proof of Lemma 6. It suffices to show that for XAi G (0, oo)"^', z = 1,..., A: and r G (0, oo) 


k 

P(Xai <xai,...,Aaj, <xaY) = <XAi)- 

i=l 

Using the notation n = YljieAi '“h = fo^ L ^ and 1) = max^.g^i Uj^Xj^, 

i = 1,..., A:, we can rewrite this equality in the form 

k 

P (Ti < rf\ ...,Yk< r-i) = [JP (U* < r"^) , 

1 = 1 


where the random vector (Ti,... , Y^) is simple max-stable (de Haan 1978) and has pair¬ 
wise independent components due to our assumptions. Hence, by Berman (1961/1962) 
Theorem 2, the 1) are jointly independent, which entails the relation above. □ 


Proof of Corollary 2. | implies Aa^^ ih Aa^^ I A; H 7^ *2 

and hence Aa^^ TL Aa^^ by Theorem 1. The hypothesis follows if we apply Lemma 6 to 
the Aa^, z = 1,..., A:. □ 
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